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Abstract 

The chaotic properties of Newton-Leipnik system are discussed from the view point of 
strange attractors. Previously, two strange attractors of this system were illustrated 
which occured from two different initial conditions under the same parameter condi- 
tion. It is found that above system also exhibits multiple attractors under different 
parameter values but same initial condition and we have shown the existence of three 
other strange attractors with varying dimensionality under different parametric con- 
ditions. The properties of these attractors are then analyzed on the basis of Lyapunov 
exponents, power spectra, recurrence analysis and peak-to-peak dynamics. The peak- 
to-peak dynamics relies on the low dimensionality of the chaotic attractor and allows 
to approximately model the system. Peak-to-peak plot along with return-time plot 
are then effectively used to solve the optimal control problem of the system which 
reverts the system to a periodic situation. 
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1 Introduction 



Beginning with the pioneering work of Lorenz on convectional instabihties [1] , there 
has been a lot of research work on systems of differential equations which yield 
bounded nonperiodic trajectories converging to an attractor of rather complicated 
nature. A chaotic attractor has a kind of "hank of string" like appearance and does 
not have a integral dimension. The chaotic trajectory is not one-dimensional as it 
does not return over itself, nor is it two-dimensional because it is not a surface uni- 
formly covered by trajectories, nor is it a volume uniformly filled by the trajectories 
to become a three-dimensional structure. 

In this communication, we have studied the Ncwton-Lcipnik chaotic system 
[2] which possesses more than one attractor. Newton and Leipnik obtained the set of 
differential equations from Euler rigid body equations which were modified with the 
addition of a linear feedback. Two strange attractors starting from different initial 
conditions and same parameter conditions were obtained. B. Marlin established the 
existence of closed orbits which were not asymptotically stable [3]. Chen et al.[4] 
devised a stable-manifold-based method for chaos control and synchronization. H. 
Richtcr studied the stabilization of a desired motion within one attractor as well as 
taking the system dynamics from one attractor to another [5]. 

We have shown that when the parametric conditions are varied keeping the 
initial conditions intact, the system showed the presence of multiple attractors. The 
chaoticity of these attractors are then analyzed through phase space analysis, power 
spcctrum[6,7] and lyapunov exponents [8,9]. Dimension of the attractor which is frac- 
tal is estimated. Owing to the low dimensionality of the attractors, the peak-to-peak 
dynamics[10,ll] is analyzed and a suitable controller [12,13]has been designed for the 
system which reduces it to a purely periodic state. 

2 Formulation 

2.1 Analysis of Chaos 

Newton Leipnik system is characterised by the following differential equations. 



This system has two strange attractors originating from the initial states (0.349, 0.0, —0.160)"^, 
(0.349,0.0,-0.180)"^ and system parameters a = 0.4,(3 = 0.175, which have been 
studied and analyzed previously [2, 4]. 

But while exploring the dynamical system with respect to the parameters we 
came across three other chaotic attractors. These three chaotic attractors emerged 
from a totally different parameter state and initial condition. The two dimensional 
projection of these attractors are shown in figure 1(a), 1(c), 1(e) and their correspond- 
ing power spectrum is illustrated in figure 1(b), 1(d), 1(f) respectively. The initial 



Xi — —ax\ -|- a;2 + 10x2X3 




X2 — —x\ — 0.4x2 + 5x1X3 

X3 = /9X3 — 5XiX2 




(2.3) 
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values of the variable, in our case, is taken as xw — 0.35, 3:20 = 0.0, x^o — 0.16. The 

parameter is kept constant at 0.17 whereas 'a' is varied. The attractor in figure 
1(a), appears for a = 0.550 whereas that of figure 1(c) is for a = 0.7689 and the last 
attractor of figure 1(e) appears at a = 0.809. Changing of the parameter 'a' causes 
a rapid change in the appearance of the strange attractor and also in its nonlinear 
properties. 

The first attractor is quite but not exactly similar to the one described in [4]. 
The second one represents much like a chaotic torus. The third one is carrying an 
entirely different structure. It is also clear from the power spectrum also that in com- 
parison to the other two the second figure shows a little bit ordered behaviour which is 
also observable in the phase space diagram. We have also studied the lyapunov spec- 
trum, which forms the most useful tool for the detection and quantification of chaos. 
Lyapunov exponents are the average exponential rates of divergence or convergence 
of nearby orbits in phase space. Any system containing atleast one positive lyapunov 
exponent is defined to be chaotic, with the magnitude of the exponent refiecting the 
time scale on which system dyanamics becomes unpredictable. The lyapunov expo- 
nents calculated for the attractors are tabulated as follows. 



TABLE I 





Ai 


A2 


A3 


Attractor 1 


0.25554 


0.03268 


-1.18911 


Attractor 2 


0.18986 


0.1039005 


-1.19466 


Attractor 3 


0.21343 


0.15888 


-1.27319 



The Lyapunov spectrum of a strange attractor is closely associated to its fractional 
dimension. Kaplan and Yorke [14-17] conjenctured that the information dimension 
df is related to the lyapunov spectrum by the equation 



df=J + 

where j is defined by the condition that 



Y^i=i Ai 
I A, 



j i+i 
Y^Xi>0 and ^ Ai < 
i=i 1=1 

The information dimension for the three attractors are obtained as, 

df, = 2.24 df^ = 2.25 d/3 = 2.29 
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It is seen that with the increase in parameter 'a', the information dimension also 
increases. Prom the lyapunov spectrum it is observed that the largest positive lya- 
punov decreases in the case of the second attractor and again rises for the third. 
So the expansion rate of the trajectories of the second attractor is less than that of 
the other two attractors whereas the third attractor displays the largest information 
dimension. 

After three attractors have been displayed with different value of parameter 'a' 
we explore the whole spectrum of the system with the variation of 'a'. This is shown 
by means of the bifurcation diagram of figure (2). There is period doubling route to 
chaos which is also depicted through phase space diagrams shown in figure (3a)- (3d). 
The first figure shows a single period with a — 0.8045. As we decrease the parameter 
value of a to 0.798 the phase space shows a double period. Further decrement of a to 
0.7944 quadrapules the period which switches to a chaotic situation after a few more 
doublings and at a value a = 0.7689. 

After this analysis on the chaotic dynamics of the system, we have taken the 
analysis of the peak-to-peak dynamics of the system which has been further used to 
design a controller for controUing the system. 

2.2 Peak-to-peak Dynamics 

It was first shown by E.N. Lorenz that for a chaotic attractor the intensity of the 
forthcoming peak and its time of occurrence can be predicted with sufficient accuracy 
by one or more curves in the plane {xi,Xi^i), where Xi denotes a maxima in the 
temporal evolution of the variable X{t). This characteristic of the attractor depends 
on the lower dimensionality of the chaotic attractor, which must be near about two. 
The subject of peak- to-peak dynamics have been discussed in the various field of 
science and engineering. Rossler system, the Chua circuit and the Duffing oscillator 
are found to display peak-to-peak dynamics under suitable parameter conditions. A 
system which display peak-to-peak dynamics can be described by a reduced order 
model which is the peak-to-peak map and a return time map. The former predicts 
the value of the next peak, while the second one gives the time interval separating a 
peak from the next peak. The methodology of obtaining a peak-to-peak plot is that, 
from the time series of a particular variable X, one collects the peak points as Xi and 
plots them Thus the plot is composed of distinct points. The attractor 

is periodic if the PPP shows only a few distinct points. The points of the PPP are 
all distinct and may show a definite curve if it is a case of torus or strange attractor 
with a dimension close to two. On the other hand for a high dimensional strange 
attractor one obtains a cloud like set. For a quasi-periodic orbit the points of PPP 
form a closed regular curve. In the similar way a return time plot can be obtained 
from the set of all pairs {xi, Ti) with = tj+i — tj. 

The peak-to-peak plot and the return time plot for the three attractors are 
depicted in figures 4(a)-4(f). The PPP and RTP for the first attractor, which occurs 
for a = 0.55 and has a fractal dimension 2.24, shows a well behaved structure as 
shown in figure 4(a) and 4(b). 4(c) and 4(d) is for the attractor with a = 0.7689 
and fractal dimension 2.25 and a distinct curve is obtained. 4(e) and 4(f) is for 
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the attractor with a — 0.809 and fractal dimension 2.29 and this is reflected in the 
PPP and RTF which has a cloud like appearance. This indicates that peak-to-peak 
dynamics is dependent on the low dimensionality of the system and is clearly visible 
from the above analysis. 

2.3 Control of Chaos 

The idea of control of chaos was enunciated by Ott, Grebogi and Yorke [22] and 
after that various methodologies and strategies have been developed theoretically 
and experimental realization of control have been achieved [23]. Control of chaos 
have been effectively observed in a magnetoelastic ribbon [24], a heart [25], a thermal 
convection loop [26], a yttrium iron garnet oscillator [27], a diode oscillator [28], 
an optical multimode chaotic solid-state laser [29], a Belousov-Zhabotinski reaction 
diffusion chemical system [30], and many other experiments. 

Peak-to-Peak plot and return time plot can be used to device optimal feedback 
controllers of chaotic systems. If the PPP is a smooth curve then one can attain a 
map from it. In our observation the attractor with a = 0.7689 showed a perfect curve 
(parabola) and it has been targeted to achieve the controlled state. This control 
procedure is adopted easily once the peak-to-peak map is derived. Our analysis has 
been considered in an ideal environment where noise is absent and the sampling 
has been made on sufficiently high frequency. In our system we introduce a control 
variable u{ij) in the third equation as follows: 

X\ — —ax\ -1-^2 + 10x2^3 (2.4) 

X2 = —xi — 0.4x2 + 5xiX3 (2.5) 

X3 = Pxs — 5xiX2 — bu(t) (2.6) 

where b is a real constant to be introduced later. For any constant control u{t) = u & 
U — [t<niin' ''^max] C R the system has a chaotic attractor. In our case for u{t) = 0.0, 
three chaotic attractor were obtained. In order to control the system, one of the 
requisites is that the control law u{t) is constant between two subsequent peaks i.e., 
u{t) = Uk V t E (tk,tk+i]- We consider the third state variable as output variable. 
A peak-to-peak map 

Xk+i = F{xk,u) 

have been derived, through simulation, for 20 equally spaced values of u in the set 
U = [0.0, 0.0095]. The reduced order model 

Xk+i = F{xk,Uk) (2.7) 

have been obtained by linear interpolation of these maps. Thus, the dynamics of 
variable y of a system can be easily described by equation (2.7). The fixed point 

X — F{x, u), 

identifies a unstable periodic orbit of the system 

x{t)=f{x{t),u{t)) 
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The output x{t) associated with the unstable periodic orbit has one peak per period 
at X. The hnear approximation of the map F around x is given by 



Xk+i - X ^ A{xk - x) + b{uk - u) 



where, 



A 



dF 
dx 



and 



dF 
du 



are two real numbers that we assume to be not equal to zero. The law for control is 
then given by 

u + h{xk — x) if I Xfe — X I < e 

Uk = < 

u ii\ Xk — X \> e 

where the scalar h must be such that \ A + hh < 1 \. So the control Uk is allowed to 
vary from {u — he) to {u + he) and 

h^-- =^ A + bh^O 
b 

is optimal because the response time of the linear approximation is minimum (i.e., 
equal to one), since, Xk+i = x. 

In the case of Newton Leipnik System, at u — 0.007 the PPM Xk+i — F{xk, u) 
has an unstable fixed point at x = 0.3 with A = —3.6 , b — 0.601. The nature of the 
controlled system with h = —A/b = 6.0 , e = 0.5 is shown in figure (5). The nature 
of the phase space is depicted in the inset that represents a single periodic orbit. 



3 Conclusion 

Three new attr actors of the Newton Leipnik system have been analyzed, which was 
found to evolve upon changing the parameter condition of the system, keeping the 
initial conditions fixed. Observation shows that they not only differ in structure and 
properties but also show varying dimensionality. Peak-to-peak dynamics have been 
obtained and the peak-to-peak plot is found to be useful in determining the controller 
function which reverts the chaotic situation to a periodic state. 
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5 Figure Captions 

Figure (1) - 2D phase space diagram and the coresponding power spectrum for 
Newton-Leipnik System, (a) & (d) : a = 0.55 (b) & (e) : a = 0.7689 (c) & (f) 
: a = 0.809. 

Figure (2) - Bifurcation diagram of the system with respect to the variation of the 
parameter a. 

Figure (3) - Period doubling route to chaos, (a) : a — 0.8045, (b) : a — 0.798, 
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(c) : q; = 0.7944, (d) : a = 0.7689. 



Figure (4) - Peak -to-peak plot and return time plot, (a) & (b): a = 0.55, (c) 
& (d): a = 0.7689, (e) & (f): a = 0.809. 

Figure (5) - Variation of X3{t) with the application of the controller showing tran- 
sition to a periodic state. Inset picture shows the phase space diagram during the 
controlled stage. 
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